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Abstract 

An exact expression for the Drude conductivity in one dimension is derived under the presence 
of an arbitrary potential. In getting the conductivity the influence of the electric field on the 
crystal potential is taken into account. This coupling leads to a systematic deformation of the 
potential and consequently to a significant modification of the charge transport. The corrections 
to the conventional Drude conductivity are determined by the configurational part of the partition 
function. The activation energy for the conductivity process is expressed by a combination of 
the free energy of the underlying equilibrium system. The electric current is calculated in the 
linear response regime by solving the Smoluchowski equation. The steady state solution differs 
significantly from the equilibrium distribution. In case of a tight binding potential the conductivity 
offers corrections depending on the amplitude of the potential. As a further application we discuss 
nanocontacts with piecewise constant potentials. The electric conductivity is corrected by the 
potential height. 
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I. INTRODUCTION 



The Drude model of electric conductivity []J offers a rather generic behavior for charge 
transport. Although the concept of the Drude conductivity seems to be old-fashioned it 
is very useful for different applications, so for chaotic transport and weak localization 2], 
magnetoresistance [3], dynamical ordering in a confined Wigner crystal J4] and optical con- 
ductivity [5[. More recently the Drude model came into the focus in describing the spin Hall 



effect [6], for a more detailed consideration see [7j. In spite of being a classical concept the 
Drude model yields the universal dc and ac conductivity and its temperature dependence 
8|. In particular, the conductivity is expressed by the charge e, the mass m and the den- 
sity n of the charge carriers as well as a relaxation time r via od = ne 2 rjm. In general 
the parameter r reflects the microscopic origin of the transport process such as different 
scattering mechanisms. In the Markov approach the relaxation time is given in terms of a 
simple damping constant 7 = r _1 . However, the Drude conductivity should be influenced 
and corrected by the underlying periodic crystal potential. And even that point we would 
like to discuss in the present paper. Here the Drude conductivity is generalized by including 
the lattice potential. The analysis is based on the Smoluchowski equation [9| which yields 
likewise the probability density for the electrons in the lattice. As the result we get an exact 
expression for the charge conductivity o = j/E in the limit E — > 0. The final relation is 
illustrated for realistic model potentials. Especially, the approach enables us to study the 
conductivity of nanocontacts by assuming piecewise constant potentials. However to find 
out the conductivity one is confronted with the problem, that the applied external electric 
field E exhibits a feedback to the potential. The potential is deformed under the influence of 
the field which leads to a modification of the motion of the charge carriers. Different to the 
quantum approach the carriers can overcome barriers or traps not by tunneling processes 
but by thermal activation. 

II. MODEL 

The classical equation of motion for the charge carriers under consideration of thermody- 
namic fluctuations are 

7 x(t) = ^ + V (t) with (r 7 (t)7 7 (t / )) = 2D5(t-t / ). (1) 
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Here 7 is a damping parameter and D is the strength of the Gaussian process. Both of them 
model the interaction with a heat bath. The Hamiltonian 7i is given by 

m 



H 



-v + U (x) — eEx . 



(2) 



where E is a homogeneous electric field and U(x) is the periodic crystal potential which 
is originated by all other particles of the system, but one can also include other potentials 
such as defect potentials. Although in the Hamiltonian the contribution of the potential 
and the electric field energy are well separated the potential may modify by the field as 



demonstrated. The Smoluchowski equation equivalent to eq. (ED) reads |9|] 



dp(x,t) dj p (x,t) 



dt 
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with j p (x,t) 
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p(x,t) . 



Inserting the Hamiltonian eq.(j2J) the probability current j p can be written as 



-D 



dp(x, t) 
dx 



1 



dUjx) 
dx 



eE 



p{x,t) . 



(3) 



(4) 



Here p(x, t) is the probability density In case the system is coupled to a heat bath with 
temperature T the equilibrium solution is given by the 



Nexp(-pH) (3~ l = k B T. 



(5) 



provided the Einstein relation D = ksT/'ym is fulfilled. M is a normalization constant. 
Notice that the equilibrium distribution satisfies j p = whereas a steady state solution 
p s (x) obeys the weaker condition d x j p = 0. It is easy to see that the equilibrium solution p e 
has a physical meaning only in case of a vanishing electric field. Otherwise, the probability 
density p e is due to the presence of the electric field not normalizable. If there exist an 
equilibrium state for nonzero field all charged particles would accumulate at x — > —00. It 
should be remarked that the same situation occurs in a quantum treatment of the problem 
where one can not find a wave function normalized in the whole space. However, different 
to the quantum case the statistical approach allows a steady state solution describing the 
physical situation adequately. The electric current is given by 

iVe 



j d = ne{v) 



dxj s (x) . 



(6) 



Here, iV is the total number of charge carriers and n is the corresponding density. From here 
we conclude immediately that the equilibrium distribution (jSJ) is an improper distribution 
function because the current is zero due to the symmetry TC(v) = TC(—v). 



III. ELECTRICAL CONDUCTIVITY 



In order to obtain the electrical conductivity one has to apply the steady state solution of 
Eq. fl3j). For a zero potential U = we find the constant steady state current by integration 
of Eq. flU 

U = ^ \d( P (0) - p{L)) + — [ L p{x)dx] . (7) 
L L l m Jo 

N is the total number of charge carriers and L is the size of the system. Taking into account 
periodic boundary conditions and the normalization condition for p(x) one obtains 
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TIC 

j e i = a D E with a D = . (8) 

rwy 

The quantity n = N/L is the particle density. The last relation is nothing else than the 
Drude conductivity where the relaxation time is given by the inverse damping parameter 7. 
For a non-zero potential the relation in Eq. (jSJ) will be modified accordingly. Directly by 
integration it follows 

Ne f L 7 dU(x) , . Ne I ' dU\ 

jei = a D E / dx— - — p s {x) = a D E ( — ) . (9) 

rwyL J dx m^L \ dx / 

The average is taken using the field dependent steady state probability distribution function 
p s (x) = p s (x;E). To find out the contribution induced by the electric field we need the 
steady state solution for the current which obeys d x j p (x) 7^ 0. Because the steady state 
solution p s depends likewise on the electric field we have to incorporate the influence of the 
electric field on the averaged crystal potential, compare Eq. 0. Since we are interested in 
the conductivity we need the steady state solution for a weak electric field E only. Therefore 
we make the ansatz 

Ps {x] E) = Af(E)e~^ [1 + E<j>{x)\ , (10) 

with an arbitrary function <p{x) and a normalization factor M{E). Inserting Eq. (1101) in 
Eq. (HD we find 

X) - 6 h ex V (U(x)/k B T). (11) 



dx k B T EDM{E) 
The current j p , defined in Eq. (J3J), disappears in for zero field. Therefore we write for small 
electric field j = aE with the unknown parameter a which is determined below. The solution 
of the Eq. (TUD is 

pt n f x 



Imposing periodic boundary conditions it follows 

eL a rL 



dx exp(U(x)/k B T) . 



(13) 



k B T DJ\f{E) Jo 

According to Eq. Q and Eq. (jTUl) as well as periodic boundary conditions the contribution 
to the conductivity is 



dU\ f L dU r 1 



.vMd<t>{x) 



e k b j 



dx 



dx . 



Using Eqs. (jlOjh pip and (jl3j) one finds 



^ ) = k B TM{E)E 



aL 



_k B TAf(0) DAf(E) 



eE 



Z L [U]Z L [-U] 



with 



Zl[U] 



dx exp([7 \x) I k B T) 



+ 0{E 2 ) (14) 



(15) 



Inserting the result in Eq. ([9]) our final relation for the conductivity reads 



L' 2 



a = a D - 



Z 2 [U = 0] 



(16) 



Z L [U]Z L [-U] " u Z L [U]Z L [-U] ' 

Because the last relation includes the ratio of the partition functions one can also insert the 
total partition function of the system without the electric field part. 

,2 



Z L [U] =Trexp[- 



mv 



+ U(x)} . 



(17) 



The final relation offers the duality property namely a symmetry against U(x) — > —U(x). 
Thus, in our approach one can not distinguish the barrier and the trapping problem. Eq. (TT6"j) 
can be also rewritten in the form 

a = a D e- EA/kBT with E A = F[U] + F[-U] - 2F[U = 0], F = -k B T\nZ . (18) 

A typical non-equilibrium quantity as the activation energy Ea is completely expressed 
by the equilibrium free energy F. Let us note that the Einstein relation is fulfilled by the 
equilibrium distribution when the system is coupled to a heat bath. Otherwise in the present 
approach we apply the steady state solution. In that case we can set D'jm = €q , where eo 
is a characteristic energy of the system, for instance the ground state energy of a quantum 
model. 
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IV. MODEL POTENTIALS 



Now let us illustrate our approach in two examples. Firstly we consider the tight binding 
potential, dotted line in Fig. [TJ 

2ir 

U(x) = Uq [1 — cos(gx) ] q= — . 

The electric field causes a systematic shift of the potential due to Eq. ( fT4l) . As the result 
the change of the probability density for the position p s (x; E) according to Eq. ( TTDT) is 
shown in Fig. [U With increasing strength of the electric field E the shift becomes more 
pronounced. The shift of the potential or the probability, respectively is so organized that a 
constant current is maintained. For a small field in terms of Uq/L the charge carriers follow 
immediately the potential. The largest probability to find the electron is at the minimum 
of the potential. When the field strength is enlarged then the probability density for the 
position is shifted considerably. For very high field the charge carriers are actually not 
influenced by the potential. In that case the probability density is nearly constant as shown 
in the last graph in Fig. [TJ The analytical calculation for the conductivity yields according 
to Eq. flUD 

' = ww- y < 19 > 

Here Io(y) is the Bessel function. In the low temperature case Uq ^> k^T the conduc- 
tivity behaves as a ~ o D exp(— 2U /kBT), i. e. the activation energy due to Eq. (TT5]) is 
dominated by the amplitude of the periodic potential Ea — 2U . In the opposite case of 
high temperatures it results Ea — Uq/21zbT. The conductivity offers a small correction 
a ~ ct.d[1 — (Po/hsT) 2 /2]. The conductivity is depicted in Fig. [2] and the activation energy 
in Fig. [3j 

A second illustration is given by a sequence of N nanocontacts modeled by a piecewise 
constant potential, U(x) = Uq < x < a and U = —Uq a < x < 2a and periodic contin- 
uation with L = Na. The number of barriers and the number of traps are both N/2. The 
conductivity follows from Eq. ( |T6l) to 

— l + col©)- (20) 
The activation energy can be expressed by 

EA = kBTI Ji_±^2lMl) (21) 



The limiting cases are 

-^4 e -W if U ^k B T; — if U <k B T. 

a D o D \ k B T J 

For an infinite high barrier the material becomes an insulator, whereas for zero potential 
the conventional Drude conductivity appears. The conductivity for the piecewise constant 
potential is shown in Fig. [2] as the dashed line. The conductivity increases continuously with 
increasing temperature. The higher the temperature in comparison to the potential height 
Uq the more charge carriers are able to overcome the barrier. For rather high temperature 
the conductivity becomes constant, i. e. the influence of the potential is negligible. The 
same situation is observed for the cos-potential (full line in Fig. El The increase of the 
conductivity is accompanied with a decrease of the activation energy Ea- The behavior of 
Ea for both the piecewise constant potential (dashed line) and the cos-potential is shown 
in Fig. [31 For practical purposes let us consider the case that there exists only one barrier 
of height U B and width A. The length of each of the two input leads is assumed to be / 
and the corresponding heights are Uq. Introducing the dimensionless ratio k = A/21 the 
conductivity can be written as 

where a m is the minimal conductivity realized for k — 1. It holds 

2od 

° m 1 + cosh(t/ B - U )/k B T) 

For k < 1 one gets 



a 

— ~ 1 — Ak 
ctd 



.Cm 

In the same manner one can find the conductivity for any other potential. 



V. CONCLUSION 



In the present paper a simple classical model for the electrical conductivity is proposed 
and solved in one dimensions. Regarding nanocontacts and nanowires such one dimensional 
models are in the focus of recent interests. Here we get an exact expression for the con- 
ductivity where systematic contributions originated by the underlying periodic potential 
are taken into account. To be specific the potential is subjected to a systematic alteration 
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caused by the applied external field. The probability distribution for the position of the 
charged particles is modified by the field. Simultaneously the charge carriers are subjected 
to an effective averaged potential which is likewise modified by the field. The alteration is so 
organized that a constant current is maintained. The contribution of the crystal potential 
is significant for low temperatures whereas in the high temperature limit the influence of 
the potential is weak and the conventional Drude conductivity is recovered. Although our 
approach is based on a classical stochastic one the conductivity can be also calculated for a 
piecewise constant potential which models nanoconducts or nanowires. The present analysis 
is restricted to the one- dimensional case. The three dimensional case is more complicated 



because the steady state solution has to satisfy the condition V • j p = w 



rich is fulfilled 



by j P = V x A with an arbitrary vector field A. In a forthcoming paper 10[ we study this 
problem in view of the spin Hall effect applying a similar approach. We are able to estimate 
the contributions of the potential to the conductivity. 

Let us finally remark that a more general expression for the conductivity is obtained when 
the Einstein relation D^ym = ksT is not longer fulfilled. Then for example the conductivity 
for the piecewise constant potential reads 

2 



cr = a D - 



1 + cosh(^a' 



whereas the characteristic energy scale eo can not be determined within our approach. How- 
ever it is appropriate to identify e with the ground state energy in presence of the potential. 
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FIG. 1: Periodic potential U(x) (dotted line) and probability distribution of the position p(x) for 
zero field (dashed lines) and with field (full line). The parameter from left top to right bottom is 
the ratio = 0.2; 1.0; 5.0; 25.0; the position x is given in units of L/2tt . 
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FIG. 2: Electric conductivity a for the cos-potential (full line) and the piecewise constant potential 
(dashed line) . 
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